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Caught in the Act: 

Discovery of a Physical Quasar Triplet 
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ABSTRACT 

We present the discovery of a triplet of quasars at z ps 1.51. The whole system is 
well accommodated within 25" (i.e., 200 kpc in projected distance). The velocity 
differences among the three objects (as measured through the broad Mg II emission 
line) are < 1000 km s -1 , suggesting that the quasars belong to the same physical 
structure. Broad band NIR images of the field do not reveal evidence of galaxies or 
galaxy clusters that could act as a gravitational lens, ruling out the possibility that 
two or all the three quasars are multiple images of a single, strongly lensed source. 
QQQ J1519+0627 is the second triplet of quasars known up to date. We estimate 
that these systems are extremely rare in terms of simple accidental superposition. 
The lack of strong galaxy overdensity suggests that this peculiar system is harboured 
in the seeds of a yet-to-be-formed massive structure. 

Based on observations collected at the La Silla Observatory with the New Tech- 
nology Telescope (NTT) of the European Southern Observatory (ESO) and at the 
Calar Alto Observatory with the 3.5m telescope of the Centro Astronmico Hispano 
Aleman (CAHA). 
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1 INTRODUCTION 

In the current hierarchical cosmological paradigm, galaxy 
mergers and interactions are a major route to galaxy 
formation. The strong tidal fields experienced during 
close encounters between gas rich systems trigger large- 
scale nuclear inflows that can ultimately start the ac- 
tivity of the supermassive black holes (SMBHs) residing 
in the central regions of the interac ting systems (e.g., 
|Pi Matteo, Springel. fc Hernquistl 120051 ) . The existence of 
multiple, simultaneously active, SMBHs represents a key ob- 
servational test for this evolutionary scenario and for our un- 
derstanding of the processes regulating the quasar (QSO) ac- 
tivity and the co-evolution of SMBHs with their host galax- 
ies. 

A number of studies demonstrated the connection be- 
tween QSOs and ultra-luminous infr ared galaxies that are 
mostly found in merging systems (e.g.. lCanalizo fc Stockton! 
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l200ll ) and confirme d the signature of re cent merger events 
in QSO hosts (e.g.. [Bennert et al.l 120081 ). A noticeable ex- 
ample is the discovery of a spatially resolved binary QSO 
(projected sepa ration pd = 21 k pc) clearly hosted by a 
galaxy merger (G reeii et al.l |2010T ). The search for candi- 
date binary systems of QSOs up to redshift ~ 4 unveiled 
an enhanced QSO clustering signature at small scale (< 
50 kpc) relative to the simple extrapolatio n of the larger 
scale two-point correlation funct i on (e.g., Hennawi et al.l 
l200rj,l2010l;lMvers et al.ll2007l,l2008l;lKavo fc Oguril2012l ; but 



Padmanabhan et al.l 120091 : IShen et all |2010| ). This can 



be interprete d in terms of interact i ons that trigger the QSO 
activity (e.g. jHennawi et al.|[200tj ; iHopkins et al.ll2008l ). Al- 
ternatively, the small scale excess could be a simple mani- 
festation of the clustering properti es of the halos that host 
QSOs (e.g., lRichardson et alj|20l5 ). 

Little is known about systems with more than two 
physically associated QSOs which are expected to be even 
more elusive objects. The only physical QSO triplet re- 
po rted so far is QQQ J1 432-0106 at z = 2.076 observed 
bv lDiorgovski et al.l ( 2007h . At smaller projected separations 
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and lower luminosities, iLiu. Shen. fc Strauss! (|201lh recently 
discovered a tri ple AGN in the galaxy S PSS J 1027+ 1749 at 
z = 0.066, and Sc hawinski et al.1 (|201 lh serendipitously ob- 
served three low mass (< 1O 7 M0) accreting black holes in 
a galaxy at z = 1.35. 

In this Paper we present the discovery of 
QQQ J1519+0627, the second triplet of QSOs that is 
known up to date. This is the first result of our ongoing 
systematic search for QSO triplets among the photometric 
and spe ctroscopic database of the Sloan Digital Sky Survey 
(SDSS; lAihara et all l201ll ). The manuscript is organised 
as follows. In Section [5] we describe the spectroscopic and 
photometric data collected in our study. The possible 
interpretations are discussed in Section [3] and we draw our 
conclusions in Section [3] 

Throughout this paper we consider a concordance cos- 
mology with Ho = 70kms _1 Mpc -1 , fi m = 0.3, and Q.a = 
0.7. All the quoted magnit ude are ex pressed in the AB stan- 
dard photometric system (|OkeHl974h . 



2 SELECTION AND FOLLOW-UP 
OBSERVATIONS 

QSO multi plets observed at clos e projected separations are 
rare (e.g., [flennawi etafl l2006h . Large spectroscopic sur- 
veys often fail to detect close QSO systems due to the 
fiber collision limits. For example in SDSS it is not pos- 
sible to obtain the spectrum for both sou rces in a pair 
with separation < 55" within a single plate (|Blanton et al.l 
2003). To overcome this limitation, we started a programme 
to search for close QSO tri plets taking advantag e of the 
large photometric sample of iRichards et al.l l|2009l ). Three 
QSOs are considered a candidate triplet if (i) at least one 
of them has spectroscopic redshift, (ii) the other two re- 
side within 500 kpc from it, and (iii) they have coinci- 
dent (within the uncertainties) photomet ric redshift (es- 
timat ed with the procedure suggested by IWeinstein et al.l 
l2004h . We selected 13 triple QSO candidates for spectro- 
scopic follow up, in the redshift range 0.3 < z < 2.2, and 
with an average projected separation of ~ 300 kpc. Fur- 
ther details on the selection procedure will be provided in 
a forthcoming paper (Farina et al., in preparation). Here 
we present the spectroscopic observations obtained for the 
first followed-up target among our candidates. In particu- 
lar, we report our discovery of two QSOs within ~ 25" from 
the spectros copic QSO SPSS J151 947.3+062753 (hereafter 
QQQ1519A, ISchneider et"afll20ld ) and located at a similar 
redshift (radial velocity difference AV < 1000 kms -1 ). 

To provide a rough estimate of the expected number of 
such systems, we derived the QSO three-point correlation 
function (£) from the amplitude of the two-point correlation 
function (£): C123M = Q [£12 (r) +63(0 + £31 M l, where 
Q « 1 fe.g.. |Peeble9ll993l ; |Pio~rgovski et alj|2007h . Assum- 
ing the va lues of the projected two-point correlation func- 
tion from iHennawi et al. |2006T l and integrat i ng th e QSO 
luminosity function presented bv lCroom et al.l |2009l ') above 
the magnitude range spanned by the photometric selected 
QSOs (i.e., from M g « -28.0 to » -25.5) we calculated that 
in a perfect survey, with no source of incompleteness, given 
a QSO, the probability of finding two companions within 
500 kpc is pabc ~ 10~ 8 . Out of the ISchneider et al l (|2010h 



catalogue, which consists in ~ 106000 spectroscopically con- 
firmed QSOs in the SPSS Pata Release 7 footprint, we thus 
expect that ~ 0.002 objects have two companions. The mere 
observation of one (or more) systems of this kind would sub- 
stantially strengthen the arg ument of small-sca le enhance- 
ment of QSO clustering (e.g.. iMvers et alj|2007h . 

In Table[T]we list the properties of the QSOs belonging 
to our newly discovered triplet, labelled QQQ J1519+0627, 
and, for comparison, of the components of the only other 
triple QSO system known so far: QQQ J1432-0106. In the 
following paragraphs we describe the procedures and the 
results of the analysis of the spectroscopic and photometric 
data collected on QQQ J1519+0627. 



2.1 Long— slit spectroscopy 

Spectra of QQQ1519B and QQQ1519C were gathered 
with the ESQ Faint Object S pectrograph and Camera 2 
(EFOSC2; iBuzzoni et all 1 19841 ) mounted on the New Tech- 
nology Telescope (NTT) in La Silla (Chile). We performed 
long-slit spectroscopy on February 17, 2012, using a slit 
width of 1" and setting the Position Angle to —49.7°, so 
that both the QSOs were observed simultaneously. Grism 
#16 was used in order to continuously cover the wavelength 
range 6015 - 10320 A with a spectral resolution A/AA « 550 
as measured on the night sky emission lines. Five frames of 
900 s were acquired, for a total exposure time of 75 min on 
source. Standard IRAlf3 tools were used in the data reduc- 
tion (bias subtraction, flat-fielding, wavelength and flux cal- 
ibration, spectra extraction). As flux calibrator we observed 
the spectrophotometric standard star Hiltner 600. Typical 
residuals in the wavelength calibration are ~ 1 A. 

For the analysis of the QSO spectra we fo llowed 
the procedure pres ented in iPecarli et al. I (I20101JI and 
|Pe Rosa et al.l |201ll ). Namely, we model the QSO spectra 
with a superposition of: (i) a power-law non-thermal com- 
ponent; (ii) the emission o f the host galaxy (ad opting the 
elliptical galaxy model by iMannucci et al.l 1200111; ( iii) the 
Fell multiplets (from the template of IVestergaard fc Wilkesl 
l200ll and from our own original spectrum of IZwOOl), and 
(iv) the broad emission lines fitted with two Gaussian c urves 
with the same peak wavelength fsee lPecarli et al . 2008). The 
two NTT spectra and the one from SPSS are shown in Fig- 
ure [1] A bright emission line, identified as Mgll, is observed 
in all the spectra at ~ 7000 A. The line identification is 
supported by the lack of other bright emission lines in the 
observed spectra, and by a tentative detection of the iron 
multiplets around the line (see our fit results in Figure [T]). 
These pieces of evidence firmly pin down the redshifts of the 
QSOs to 1.504 ± 0.001, 1.513 ± 0.003 and 1.506 ± 0.003 for 
A, B and C respectively, where uncertainties are given by 
the positional accuracy of the Mgll line centroids. 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 1. Properties of the triple QSO systems known to date. For each QSO, we list: identification label (id), position (RA, DEC), 
redshift (redshift), magnitude of the components (z, J, and H), angular ( Afl) and projected (pd) s eparations at the redshift of the system, 
and radial velocity difference (AV). Data for QQQ J1432-0106 are from lDiorgovski et al.l l|2007T ) and the SDSS database. 
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redshift 


z 
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QQQ J1519+0627: 


























- QQQ1519A 


15:19:47.3 


+06:27:53 


1.504+0.001 


18.86 


18.81 


18.55 


Afl(A-B)= 


23.5 


pd(A-B)= 


198 


AV(A-B)= 


1100 


- QQQ1519B 


15:19:45.7 


+06:27:52 


1.513+0.003 


21.23 


20.97 


20.25 


A0(B-C)= 


3.7 


pd(B-C)= 


31 


AV(B-C)= 


850 


- QQQ1519C 


15:19:45.9 


+06:27:49 


1.506+0.003 


21.20 


20.69 


20.13 


A0(C-A)= 


21.1 


pd(C-A)= 


178 


AV(C-A)= 


250 


QQQ J1432-0106: 


























- QQQ1432A 


14:32:29.2 


-01:06:16 


2.076 


17.15 


16.32 


15.88 


A0(A-B)= 


5.1 


pd(A-B)= 


42 


AV(A-B)= 


280 


- QQQ1432B 


14:32:28.9 


-01:06:13 


2.076 


20.42 


20.00 


19.27 


A0(B-C)= 


3.6 


pd(B-C)= 


30 


AV(B-C) = 


100 


- QQQ1432C 


14:32:29.2 


-01:06:12 


~2.08 




21.66 


21.17 


A0(C-A)= 


4.3 


pd(C-A)= 


36 
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Figure 1. Spectra of the three QSOs of QQQ J1519+0627 (grey), binned by 2 pixels. The results of the fitting procedure on Mgll and 
Fell lines are plotted in blue (see text for details). The positions of the Mgll peaks are marked with orange dashed lines. 

2.2 Broad band photometry 

We gathered NIR broad band ima ges of the 
QQQ J1519+0627 field using Omega2000 (|Kovacs et aD 
|2004| ) at the 3.5m telescope in Calar Alto (Spain). Ob- 
servations were performed on March 8, 2012 as Director 
Discretional Time. We observed a region of 15' x 15' 
around the system in z, J, and H bands, i.e., sampling 
the rest-frame u, g, and r at z = 1.51. We adapted usual 
jittering observing strategies in order to effectively subtract 
night sky emission. We collected 180 x 15 s frames in z 
(total integration time: 45min), 18 x 100 s frames in J 



(total integration time: 30min), and 36 x 50 s frames in H 
(total integration time: 30min). Images were processed with 
our package of NIR image processing based on standard 
IRAF tasks. We compute t he astrometric solu tion using the 
astrometry.net software ijLang et al . 2010). Photometric 
calibration is achieved by comparing the photometry of field 
stars observed in our images with the fluxes reported in the 
SDSS (z-band) and 2MASS (J- and H-band) catalogues, 
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using the following filter transformations: 
zrm = zsdss — 0.05 (i — z) SDSS 

Jf!2fc = J2MASS.AB + 0.11 (J — H) 2MAggAB 

H n2fe = H 2M ass,ab - 0.02 (J - H) 



(1) 

(2) 
(3) 



2MASS,AB 

Equations Q]-[3] are derived by computing the spectral mag- 
nitudes of template O- to M-type stars. 

The seeing during the observations was 2" (as measured 
on the final images), and the 5<r limit magnitudes (com- 
puted from the rms of sky counts over a seeing area) are 
23.90, 23.38 and 22.25 in z, J and H respect i vely. 

We use SExtractor (1 Bert in fe Arnoutslll99tjt ) in order 
to identify and catalogue the sources in our images. Within 
the limit fluxes of our observations, we are able to study the 
environment of the QSO system down to M*(z) + 1, where 
M* (z) is the rest-frame B-band characteristic luminosity of 
galaxies at redshift z as derived from the fit of the galaxy 
luminosity function |llbert et al.l 12005), In order to distan- 
gle between Galactic and extragalactic objects we rely upon 
different colour-colour diagrams, also including sources from 
the SDSS photometric database (see Figure [2]). 

We classified a source as star if it is not resolved 
(i.e., FWHM < seeing in all the considered images), 
and it lies within the locus of main sequence stars esti- 
mated through the prediction of the Trilegal softwarfl 
jGirardi et aij|2005h . and from the distribution of the SDSS 
spectroscopically con firmed stars detected also in 2MASS 
l|Skrutskie et alj|2006h . Our colour cuts are a s implified ver- 
sion of those adopted bv lRichards et all |2002l ) on the SDSS 
photometric database. We also check the consistency of the 
sources not classified as stars with the reference colours of 
Elliptica l, Sc galaxies and Q SOs obt ained from the tem - 
plates of lMannucci et al.1 l|200ll ) and of lFrancis et alj (|l99lf ). 
In summary in the 15' x 15' region explored we detect 
2048 sources, 807 of which are classified as star on the ba- 
sis of their colours. This yields a number density of ~ 3.5 
stars per arcmin 2 that roughly correspond to the predic- 
tion of the Trilegal software in this region of the sky 
(~ 3.1 arcmin -2 ). The average surface density of extragalac- 
tic sources is 3.7 ± 0.2 arcmin -2 , that is consistent (within 
the uncertainties) with the 4.0 ± 0.2 arcmin -2 observed in 
the The Great Observ atories Origins Deep Survey (GOODS, 
iGiavalisco et alj|2004l ) once we consider our z-band sensitiv- 
ity limits. 



3 DISCUSSION 

In this section we discuss the possible scenarios for this pe- 
culiar system. 



3.1 Gravitational lens 

To test the lens hypothesis for QQQ J1519+0627 we con- 
sider two different configurations: in case (i) QQQ1519B and 
QQQ1519C are images of the same, gravitationally lensed 
QSO while QQQ1519A is a different QSO at similar redshift. 
In case (ii), also QQQ1519A is an image of the same QSO. 
We note that small differences are reported in the colours 
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Figure 3. Expected Einstein ring radius for a lensing galaxy 
with various (observed) z-band magnitudes, as a function of 
redshift. The horizontal line shows the half-separation be- 
tween QQQ1519B an d C. The elliptical galaxy template by 
Mannu cci et al . (2001) is used in order to compute k— and filter- 
corrections. We assume that the potential well of the lensing 
galaxy is well described by a SIS model with a = cr*, where 
<7* is the stellar velocity dispersion as set by the Faber-Jackson 
relation (see text for details). In the lensing scenario, a galaxy of 
zRil9.6 would be required in order to justify the observed sepa- 
ration between QQQ1519B and C. Thanks to the depth of our 
observations (zi; m PS 23.9), we can rule out the strong-lensing 
scenario for these two sources. 



and in the spectra of the three QSOs (e.g., the Mgll lines 
peak at slightly different redshifts). This already disfavours 
the lensing scenario. However, due to the limited signal-to- 
noise ratio of the available spectra, we will focus our analysis 
mostly on the images. 

In case (i) we model the potential well of the lensing 
object with a singular isothermal sphere (SIS). Multiple im- 
ages of similar brightness are obtained only along the Ein- 
stein ring, the radius of which is parametrised as 9e (e.g., 
Naravan fc Schneider! Il99d : IChieregato. Miranda. fc Jetzerl 
2007ft : 



: 4"/T 



(<t\ 2 Dls 
\c) D s 



(4) 



Here a is the velocity dispersion of the SIS, c is the speed 
of light, Dls and Ds are the angular diameter distances 
between the source and the lens, and between the source and 
the observer, respectively. The minimum value of 8e allowed 
in this system corresponds to half of the separation in the 
sky between QQQ1519B and C, i.e., l'.'85. This implies that, 
for a given redshift of the lens, we can infer a lower limit for 
a. If we adopt the working assumption that a = a* (the 
velocity dispersion of st ars in the galaxy), we can use the 
Faber-Jackson relation (Fab er fc Jackson! Il976t ) to convert 
the lower limit on a into a luminosity limit for the lens 
galaxy. We use the updated version of the Faber-Jackson 
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Figure 2. Colour-colour diagrams of the sources detected in the explored images (grey points). Objects within 60" from the 
QQQ J1519+0627 baricenter are highlighted in red. The blue triangles indicate the colours of the three QSOs of the system. Note 
that the number of detected sources changes, according to the different depth of each band. For compar ison we show the expecte d loci 
of Ellipticals and Sc galaxies at z = 1 and of QSOs at z = 1.5 from the templ ate oflMan nucci et al.l ^200ll ) and lFrancis et af | Jl99lh . The 
light blue points are the star colours simulated with the Trilegal software iGirardi*^t**aLll2005l) . and dark green contours indicate the 
distribution of the spcctroscopically confirmed stars detected in the SDSS and 2MASS in a region of 5° X 5° around the triplet, and in 
the GOODS survey. Black dashed lines show the colour cuts we imposed to separate stars and extragalactic objects. The crosses in the 
bottom-right of the figures represent typical error bars on the colours. 



relation reported by iNigoche-Netro et al.l ((2010) , 



log 



kms~ 



= (-1.208 ± 0.205) - (0.157 ± 0.009) M r (5) 



where M r is the rest-frame r-band absolute magnitude. We 
convert the observed z-band into M r using the filter and 
fc-correct ions computed b a sing o n the elliptical galaxy tem- 
plate by iMannucci et al.l (|200ll ). Figure [3] shows that a 
galaxy brighter than z=19.6 is necessary to explain the ob- 
served geometry of the system in terms of strong gravita- 
tional lensing. The presence of such a galaxy is ruled out by 
our Omega2000 images. 

Case (ii) is more extreme since wide separation lens 
systems are rare. In their search for le nsed QSOs w i th an - 
gular separation less than A9 < 20", llnada et all (|2008l ) 
find in the SDSS only 3 out o f 22 systems with Ad > 3". 
iHennawi, Dalai, fc Bode] (|2007f) estimated that in the SDSS 
QSO sample there should be only a few systems (< 4) multi- 
ply lensed by galaxy clusters with separations > 20". More- 
over, triple imaged Q SOs are not commonly observed (e.g., 
lLawrence et al.lll9&4r ). 

The three largest separation QSO lenses k nown so far 
22" 5, z 
15'.' 1 



SDSS J1029+2623 (A<9 



2006), SDSS J2222+2745 (AO 



2.197, 
= 2.82. 



Inada et al 



Dahleetal 



2012) , and SPSS J1004+4112 (AO = 14'' 6, z = 1.734, 



Inada et al.ll2003l ). have separation similar to the ones re- 



ported here. In all these cases massive galaxy clusters act 



as le ns l|Oguri et all |2004 llnada et all l200rj ; iDahle et all 
l2012h . As a zeroth order test, we estimate the velocity 
dispersion required to explain the separation of the im- 
ages assuming again a SIS profile for the lens. As in case 
(i) we consider a limit for the Einstein radius 9e ~ 10" 
that corresponds to a ~ 600 km s -1 at z m and a w 
1200kms _1 at z ~ 1. In our image the presence of such 
clusters should be apparent, but in fact the surface density 
of galaxies within 60" from QQQ J1519+0627 (in z-band 
4.1 ± l.larcmin -2 ) is consistent with those detected in the 
whole frame (3.7±0.2 arcmin -2 , see i]2.2l and Figure 2J. For 
instance, we consider the case that in our field was present 
XMMU J100750. 5+125818: a galaxy cluster located at red- 
shift z ~ 1.08 with a line-of-sigh t velocity dispersion of 
~ 600kms _1 (|Schwope et aU feoiO). If the cluster was cen- 
tred on QQQ J1519+0627 in the first 60" there should be 
at least 10 galaxies in excess to the background. This is 
in contrast to what observed in the z-band where we de- 
tect only 12 galaxies instead of 22. We note that this esti- 
mate is conservative since we considered as reference only 
the (few) spectroscopically confirmed cluster members of 
XMMU J100750.5+125818, while the cluster should be more 
populated. 

From the analysis of photometric data, we thus conclude 
that QQQ J1519+0627 is not a case of gravitation lens. 
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Figure 4. Image of the 60" around the triplet QQQ J1519+0627 
as imaged in z— band at the Calar Alto Observatory. Blue dia- 
monds are the three QSOs. The 13 sources identified as galaxies 
are shown as red squares, while the 9 stars as green circles. The 
surface densities of Galactic and extragalactic sources in this area 
is consistent with the detections in the rest of the image. 



we search for galaxies at redshift ~ 1.5 in the proximity 
of the system. In Figure [5] we show the observed colour- 
magnitude diagram of the sources in our field. We compare 
the photometry of these sources with the apparent mag- 
nitude and colours of a M*(z) galaxy, where M*(z) is the 
characteristic luminosity (in absolute magnitudes) of galax- 
ies at a given redshift z in the rest-frame B-band, as com- 
puted bv lllbert et all |2005l ). For filter and fc-correction, we 
refer to the elliptical galaxy template by iMannucci et al.l 
|200ll ). We do not find evidence of a strong overdensity of 
sources with colours consistent with those of a red sequence 
at z « 1.5. This suggests that QQQ J1519+0627 is not lo- 
cated in a rich galaxy cluster or that the red sequence is not 
yet formed. We also estimate the photometric redshifts of 
the sources in the field of QQQ J1 519+0627 with Hyper^ 
jBolzonella. Miralles. fc Pellol [200ch relying upon different 
galaxy templates (i.e., starburst, elliptical, lenticular, dif- 
ferent type of spirals, an d irregular ga l axies) and assum- 
ing the extinction law bv lCalzetti et all (|2000h . It is worth 
noting that this choice of dust extinction model has only 
marginal effects on our results. This analysis confirms that 
none of the objects detected within 60" (that corresponds to 
~ 500 kpc at the triplet redshift) and 120" (~ 1 Mpc) have 
colours consistent with a galaxy at z ~ 1.5 (see Figure [5]). 
We thus propose that the physical structure, where the three 
QSOs reside, may be caught in the act of formation (sim- 
ilarly to what reported in e.g.: iDecarli. Reynolds, fe Dottil 
|2009| ; IDecarli et ai1l2010al ; iFarina. Falomo. fc Trevesll201ll ). 



3.2 Chance superposition 

QQQ J1519+0627 is the second triplet of QSOs at simi- 
lar re dshift known, after QQQ J1432-0106 |Diorgovski et all 
2007). In Section [2] we have already shown that these are 
extremely rare objects. Here we estimate the expected num- 
ber of QSO triplets in our sample that are originated by 
chance superposition of otherwise (physically) disconnected 
QSOs. We start from the sample s of sp ectroscopically con- 
firmed QSOs by ISchneider et alj {2010) an d the sample of 
photo metrically-selected QSO candidates bv lRichards et al] 
(2009) (both based on SDSS data) and we rely upon a pro- 
cedure similar to the redshift perm utation method (e.g., 
IOsmer|[l98ll ; IZhdanov fc Surdeill200ll 'l. We assume that the 
QSOs have the same position in the sky as in the original 
samples but a random redshift is assigned through a Monte 
Carlo simulation that takes into account the redshift dis- 
tribution of our complete sample (~ 940000 objects). We 
find that ~ 0.05 QSO triplets with AV < 2000 km s" 1 and 
projected separations < 200 kpc are expected. This result is 
an upper limit for the number of chance triple QSOs since 
in our computation most of the correlation was destroyed, 
but we kept the angular distribution of the QSOs as in the 
original sample. 

3.3 A common physical structure 

In the hypothesis that the system is virialised, the velocity 
differences among QQQ1519A, QQQ1519B and QQQ1519C 
(estimated from the broad Mg II emission line redshifts) im- 
ply a dynamical mass of ~ 10 13 M©, i.e. the system would 
lie inside a massive structure. In order to test this scenario, 



4 CONCLUSIONS 

We report the discovery of a triplet of QSOs with angular 
separations between 3" 7 and 23" 5 at a similar redshift of 
z « 1.5. The three sources have relative velocities within 
~ lOOOkms -1 as derived from our measurements of the 
Mg II broad emission line and of the surrounding Fe II mul- 
tiplets in the observed optical spectra. From the analysis of 
u, g, r, i, z, J, and H broad band images we conclude that 
this system consist of three distinct sources rather than be- 
ing a multiply imaged lensed QSO. This would be the sec- 
ond triplet of QSOs disc overed so far after QQQ J1432-0106 
(|Diorgovski et al] 120071 ). We estimate that the probability 
that these systems are due to chance superposition are negli- 
gible. No clear galaxy overdensity is reported at photometric 
redshift consistent with the one of the triplet. We therefore 
propose that these systems are part of an ongoing common 
physical structure formation. 

The projected separation of QQQ1519A from the other 
two members of the system is greater than the typi- 
cal distances between interacting systems (~ 50 kpc; e.g ., 
iHennawi et al. I l2006l ; iForeman. Volonteri. fc Dottil [2009). 
This implies that the observed nuclear activity of 
QQQ1519A was probably not triggered during an inter- 
action involving the three systems at the same time. On 
the other hand, due to the smaller separation observed be- 
tween QQQ1519B and C, the triggering of nuclear activity 
by galactic interactions may be a viable scenario for these 
QSOs. More conclusive constraints on the dynamics of this 
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Figure 5. Observed colour— magnitude diagram for the objects 
detected by SExtractor in the field of QQQ J1519+0627(green 
circles). Sources selected as galaxies (see <12.2[l are shown in or- 
ange, and those that lie within 60" from the triplet are highlighted 
with red squares. The vertical dotted line marks the 5 a detection 
limit (note that SExtractor is on average more conservative than 
the statistical limit magnitude). Typical error bars are shown as 
a cross in the bottom-right of the figure. The dashed line shows 
the l ocation of an M* ( z) galaxy at various rcdshifts, as derived 
from lllbert et al. 1 J2005h. and as suming the elliptical galaxy tem- 
plate byjM^rmua^^^I ] J200 J l for fc-corrections. For reference, 
we show with light violet crosses the location of the galaxies of 
the galaxy cluster XMM XCS J2215.9-1738 (z = 1.46, velocity dis- 
persi on <7„ = 580 km s- 1 , [Stanford et al.ll2006l;lHilton et aT1l2007l. 
2009). Dark violet crosses highlight the spectroscopic confirmed 
associations. We do not find any evidence of an overdensity of 
galaxies with colours consistent with a red sequence at 2 = 1.51, 
thus suggesting that the environment of QQQ J1519+0627 is not 
a rich cluster. 



peculiar system and on the properties of its environment will 
be obtained from future higher signal-to-noise ratio spectro- 
scopic data in the NIR band and deeper photometric images. 
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